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Previous studies have demonstrated that the use of counterflowing jets can greatly reduce the 
drag and heat loads on blunt-body geometries, especially when the long penetration mode jet 
condition can be established. Previously, the authors had done some preliminary numerical studies 
to determine the ability to establish long penetration mode jets on a typical Mach 1.6 slender 
configuration, and study its impact on the boom signature. The results indicated that a jet with a 
longer penetration length was required to achieve any impact on the boom signature of a typical 
Mach 1.6 slender configuration. This paper focuses on an in-depth parametric study, done using the 
space-time conservation element solution element Navier-Stokes flow solver, for investigating the 
effect of various counterflowing jet conditions/configurations on two supersonic slender-body 
models (cone-cylinder and quartic body of revolution). The study is aimed at gaining a better 
understanding of the relationship between the shock penetration length and reduction of drag and 
boom signature for these two supersonic slender-body configurations. Different jet flow rates, 
Mach numbers, nozzle jet exit diameters and jet-to-base diameter ratios were examined. The 
results show the characteristics of a short-to-long-to-short penetration-mode pattern with the 
increase of jet mass flow rates, observed across various counterflowing jet nozzle configurations. 
Though the optimal shock penetration length for potential boom-signature mitigation is tied to the 
long penetration mode, it often results in a very unsteady flow and leads to large oscillations of 
surface pressure and drag. Furthermore, depending on the geometry of the slender body, longer jet 
penetration did not always result in maximum drag reduction. For the quartic geometry, the 
maximum drag reduction corresponds well to the longest shock penetration length, while this was 
not the case for the cone-cylinder — as the geometry was already optimized for drag. Numerical 
results and assessments obtained from this parametric study along with the recommendation for 
future implementation of counterflowing jets as a means for drag and noise reduction are detailed 
in this paper. 
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Nomenclature 

Jet exit diameter, mm (or inch) 

Model diameter, mm (or inch) 

distance away from the centerline of the body, cm 

Length, cm 

Jet penetration length, cm (or inch) 
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= Counterflowing jet nozzle Mach number 
= Freestream Mach number 
= Freestream Mach number 
= Pressure, atm 

= Freestream static pressure, atm 
= Total pressure of counterflowing jet, atm 
= Total freestream pressure, atm 
= Total freestream temperature, K 
= Cartesian coordinates 


I. Background and Introduction 

Counterflowing jets emanating from blunt bodies or bodies of revolution against a supersonic/hypersonic 
freestream have been the subject of several experimental studies, 1 ' 8 as well as a small number of computational studies. 9 ' 14 
Two different modes of jet interaction with the oncoming freestream exist, namely the short penetration mode (SPM) and 
long penetration mode (LPM), pictorially represented in Fig. 1. Both these modes have potential benefits for drag or thermal 
load reduction. Depending on the location of the jet relative to the face of the vehicle forebody, both thermal load reduction 
and drag reduction/augmentation on the body are possible. The LPM jet, when issuing from a centrally located single 
nozzle, often penetrates into the bow shock in front of the body, weakening its strength or dispersing it, creating a complex 
flow interaction. As a result, it significantly alters the flowfield around the body, thereby having a potential to be applied as 
a means to reduce sonic boom for vehicles travelling at supersonic speeds. 

Sonic boom mitigation remains an active field of research since the 1960s, 15 although the theory of sonic boom 16 ' 18 
and the methods for its prediction 19 ' 21 have been well established. In 1969, Miller and Carlson 22 explored an idea similar to 
that of counterflowing jets for sonic boom mitigation, where through the addition of heat/forcefields, they attempted to 
create a long, well-shaped phantom body that envelops the aircraft. The conclusion of their study was that such an idea 
would necessitate power requirements that are at least twice of what is needed to sustain steady flight, even if the system 
required for achieving heat/mass addition was assumed to be weightless. More recently, Fomin et al. 23 have also explored 
the concept of using counterflowing jets as a boom mitigation technique, based on their prior experience with 
counterflowing jets. 3 Their experimental study focused predominantly on SPM jets, giving limited attention to LPM jets. 
The conclusion from their study was that counterflowing jets did not yield significant boom signature reduction for the 
cone-cylinder type geometry they had investigated. 




(a) SPM: features of the stable flowfield for a central nozzle (b) LPM: features of an unstable flowfield for a central 
configuration. nozzle configuration 

Figure 1. Characterization of the flowfield structure for a jet emanating from a central nozzle against a supersonic 
freestream. 

Despite the conclusions of Miller and Carlson 22 and Fomin et al., 23 we intended to explore the potential of LPM jets 
and their application in boom and drag mitigation in a more thorough manner, using computational studies. Our improved 
understanding of the LPM jets through earlier investigations 9, 14 drove this study. In this regard, the authors of this paper had 
previously conducted a pilot computational study 24 on the use of the LPM of a counterflowing supersonic jet for sonic boom 
mitigation. In that study, the jets emanated from the nose of two slender-body geometries 25, 26 (a cone-cylinder and a quartic 
body of revolution), which were of interest to the High Speed project under NASA’s Fundamental Aeronautics Program 
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(FAP). Based on the understanding that LPM jets exist under a narrow range of flow conditions, the main focus of that 
study was to investigate the feasibility of obtaining an LPM jet issuing from these two slender-body configurations against 
low supersonic freestream conditions. Numerical computations were carried out to study a limited range of counterflowing 
jet conditions and nozzle geometries for establishing LPM jets. The results indicated that the jet penetration lengths (i.e., 
shock stand-off distance) achieved might not be long enough to gain substantial reduction in boom signature. In this study, 
we extend the previous preliminary study with a wider range of key parameters associated with counterflowing jet 
conditions, and explore the optimal condition for the longest jet penetration length. The objectives of the present numerical 
study are (i) to assess the ability to establish counterflowing LPM jet and its associated penetration length as well as 
interaction with the shock system around the body of a cone-cylinder 25 and a quartic geometry 26 at a freestream Mach 
number of 1.6; (ii) to determine the effects of key geometric parameters and flow conditions on the jet penetration length 
and drag change; (iii) to understand the correlation between the jet penetration length and drag change; and (iv) to assess the 
influence of the LPM jet on the farfield pressure signature. The numerical method used in this study is the space-time 
conservation element, solution element (CESE) method, 27-29 a time-accurate numerical framework for unstructured meshes 
designed to enforce strong flux conservation and capture shocks without a need for any ad hoc numerical tuning. The 
ability of the CESE method to correctly predict the occurrence of LPM jets, as well as its behavior, has been demonstrated 
in our previous studies. 9, 14 


II. Test Geometries and Flow Conditions 

Because of the interest of these geometries to the High Speed project, all the numerical computations of this study 
have been carried out on an 21.844 cm (8.6 in) long 6.48° (full-angle) cone-cylinder and a 5.08 cm (2 in) long quartic 
geometry, described in Refs. 25 and 26, respectively. The freestream condition tested in this study was a Mach 1.6 flow at 
an altitude of 13,700 m and 0° angle of attack, as these conditions were identified by the High Speed project to be optimal 
conditions for a potential supersonic flight. This corresponds to a static pressure of 0.1462 atm, a static temperature of 
216.65 K, and a unit Reynolds number of 7.912xl0 6 m" 1 . To facilitate injecting a counterflowing jet from the body’s nose, 
some modifications were performed to the cone-cylinder geometry, while no major modification was needed for the quartic 
geometry, as seen in Fig. 2. The cone-cylinder was truncated with different lengths at the nose to accommodate various 
nozzle exit diameters from where the jet exhausts. Three bell nozzle configurations for achieving different jet exit Mach 
numbers (2.94, 3.52, and 4.02) were used in this study to assess the jet Mach number effect, while each nozzle contour was 
scaled to fit the desired jet exit diameter. The nozzle inlet was assumed to be connected to a reservoir chamber at stagnation 
conditions of 300 K and a pressure, which yields the desired jet exit pressure and its corresponding mass flow rate. 



(a) Truncated cone-cylinder with nozzle (b) Quartic geometry with nozzle 


Figure 2. Modified geometries with nozzle included. 

III. Computational Details 

As the focus of these studies was to gain a better understanding about the ability to establish LPM jets in slender 
bodies at low supersonic freestream conditions and to determine the factors influential in obtaining longer LPM jets, the 
computational domain, shown in Fig. 3, encompasses only the near-field region. This is different from the approach adopted 
in many other sonic boom prediction studies, 30 where the computational domain extends much further into the free stream 
and the computations are mainly done using inviscid flow assumptions for simple axisymmetric geometries. The use of 
smaller computational domains permits one to more cost-effectively conduct the multiple computational fluid dynamics 
(CFD) simulations needed in a parametric study. Moreover, viscous effects are very important in this study because of the 
LPM jet being the primary focus. Through previous studies on the counterflowing jets, 9, u " 14, 24 we have understood that 
capturing LPM jets requires an extremely fine mesh. As a result, we devote one region of the mesh towards capturing the jet 
development, denoted as the jet flow region in Fig. 3, where the cells remain almost uniform in cell size, with y + (non- 
dimensional grid spacing normal to the wall, calculated based on the jet Reynolds number) varying from 0.5 near the wall, 
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to a maximum of 4.0. This region extends several nozzle exit-diameter lengths away from the nose of the body. Away from 
this region, the mesh still needs to be fine enough to predict the propagation of the pressure disturbances arising from the 
jet, but not as fine as the jet flow region. Complete details about the grid resolution requirements, and results from grid 
convergence studies were given in our previous work. 24 

The CESE-based Navier-Stokes code, ez4d , 29 developed at the NASA Langley Research Center, was employed to 
perform all the computations in this study, with the assumption of axisymmetric laminar flows. Unstructured meshes were 
utilized for all numerical computations carried out in this study. The boundary conditions employed are shown in Fig. 3. In 
particular, the time-accurate local time stepping scheme (TALTS) 31 was employed to minimize the numerical dissipation 
induced by the large disparity in grid sizes. In general, much smaller meshes are used in the region with large flow 
gradients, such as shock, boundary and shear layers, compared to the freestream region. The time step size in the TALTS 
calculations varies with the mesh size to achieve a relatively uniform Courant-Friedrichs-Lewy (CFL) number (around 0.8) 
throughout the computational domain so that numerical dissipation is minimized and the structure of jet expansion, the 
weak bow shock, and the unsteadiness of shock locations can be accurately captured. The jet Reynolds number (based on 
the nozzle exit diameter and the nozzle flow conditions) for all of the cases computed was in the range of 10,000 - 30,000. 
Per the experimental study by McLaughlin et al., 32 this falls within the low-to-moderate Reynolds number range where the 
flow was found to be in a laminar or at most transitional regime. Hence, the laminar flow assumption made here is 
reasonably justified. 



Figure 3. Computational domain and boundary condition details. 

IV. Numerical Results from the Parametric Study 

Some key parameters, indicated by earlier studies on counterflowing jets, 1 ' 3, 24, 33 were explored towards achieving 
LPM jets (issuing from slender-bodies under low supersonic freestream conditions) that are capable of pushing the leading 
shock as far from the body as possible. For the cone-cylinder geometry, these parameters include (1) jet/nozzle exit Mach 
number (Af.), (2) jet/nozzle exit diameter (J), (3) jet/nozzle total pressure (P ) or nozzle-to-freestream pressure ratio (NPR), 
and (4) jet-to-base diameter ratio (d. IdJ. With the jet mass flow rate being directly tied to the nozzle total pressure, its 
variation is equivalent to the variation of nozzle total pressure or NPR. For the quartic geometry, the effect of d. is the same 
as d.ld m since d m (= 5.733 mm, or 0.2257 in) is fixed (see Fig. 2b). Otherwise, parameters (l)-(3) mentioned above were 
explored for the quartic geometry also. In our previous study, 24 the effect of nozzle divergence angle on jet development 
was found to be very small for semi -divergent angles up to 20°, and thus, was not investigated in this study. The results, 
regarding the drag change and the jet penetration lengths, are presented for the cone-cylinder followed by those of the 
quartic geometry 

A. Cone-Cylinder Geometry 

With the total drag (pressure and viscous drag) being used as one of parameters for assessing the effect of the 
counterflowing jet, the drag of the baseline cone-cylinder was first computed. The calculated drag for the section of the 
cone-cylinder simulated converged to 1.328 N and this value was employed to calculate the percentage of drag change for 
the use of the counterflowing jet. A summary of the various counterflowing jet configurations studied is given in Table 1. 

a. Effect of Jet/Nozzle Exit Mach number 

To examine the effect of jet/nozzle exit Mach number (Af.), three different nozzle contours (i.e., throat-to-exit area 
ratios) were used with the same nozzle exit diameter (d. = 0.635 mm, or 0.025 in) and jet-to-base diameter ratio ( d./d m = 
1/4), while the nozzle total pressure is varied to obtain a trend. Three bell-shaped nozzles were selected to provide jet exit 
Mach numbers of 2.94, 3.52 and 4.02, respectively. The corresponding range of nozzle total pressures is different for each 
of the nozzles, so that LPM can be established and the produced jet mass flow rates are comparable in all cases. The results 
indicate that a longer jet penetration length (L p ) can be achieved with higher Af. and NPRs (thus, the jet flow rates), as 
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shown in Figs. 4-6. The error bars shown in these plots and many of the other plots included in this work represent the 
maximum and minimum values of the shock standoff distance and drag change, observed over the duration of the 
simulation, due to the inherent unsteadiness of the LPM jet. 


Table 1. Summary of counterflowing jet configurations studied for the cone-cylinder geometry 


d. 

j 

d.ld 

j m 

M 

j 

NPR 

f’tj (atm) 

0.3175 mm (0.0125 in) 

1/4 

2.94 

34.2-88.9 

5-13 

0.635 mm (0.025 in) 

1/4 

2.94 

34.2- 102.6 

5-15 

0.635 mm (0.025 in) 

1/4 

3.52 

75.2-205.2 

11-30 

0.635 mm (0.025 in) 

1/4 

4.02 

116.8-321.5 

17-47 

0.635 mm (0.025 in) 

1/8 

2.94 

34.2- 102.6 

5 -15 

0.635 mm (0.025 in) 

1/16 

2.94 

34.2- 102.6 

5-15 

1.27 mm (0.05 in) 

1/4 

2.94 

20.5-88.9 

3-13 

1.27 mm (0.05 in) 

1 

4.02 

116.3-403.6 

17-59 



Figure 4. Effect of jet exit Mach number on jet penetration length with various NPR for the cone-cylinder geometry 
(i d . / d m = 1/4, and d. = 0.635 mm or 0.025 in). The error bars correspond to maximum and minimum values of a 
variable observed over the duration of the simulation. 



Figure 5. Effect of jet exit Mach number on jet penetration length and drag reduction with various jet total 
pressures for the cone-cylinder geometry ( dj !d m = 1/4, dj = 0.635 mm or 0.025 in). 
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Jet flow rate 10 4 (kg/s) Jet flow rate 10 4 (kg/s) 


Figure 6. Effect of jet exit Mach umber on jet penetration length and drag reduction with various jet flow rates for 
the cone-cylinder geometry (d. Id m = 1/4, and d. = 0.635 mm or 0.025 in). 
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M: 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 4 4.2 4.4 



M. = 4.02 
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Figure 7. Instantaneous Mach contours showing the effect of jet exit Mach number on the jet development for the 
cone-cylinder geometry. The x- and y-coordinates are given in inches. The nose tip of the unmodified cone-cylinder 
geometry was at x = 0. 
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The Mach number contours, corresponding to the longest L p for each M. for a given nozzle total pressure is shown in 
Fig. 7. As can be seen from Fig. 7, even under conditions with the longest L , the jet penetration is unable to compensate for 
the truncation in the length of the cone-cylinder (in Fig. 7, the tip of the unmodified cone-cylinder was at x = 0), which was 
necessary for implementation of the counterflowing jet nozzle. Furthermore, for a given M., L increases over a range of 
NPRs (and jet flow rates) to an optimal value and then decreases — an indication of the jet switching from LPM to SPM. The 
data points corresponding to NPR = 102.6 for M. = 2.94, and NPR = 205.2 for M. = 3.52 in Fig. 4 are indicative of the 
transition to SPM. The existence of a range of pressures at which the LPM jet regime can be established offers some 
flexibility, in case the freestream condition changes continuously. This aspect could become important, especially during 
vehicle acceleration/deceleration at supersonic speeds. In addition, drag reduction increases monotonically with NPR (and 
jet flow rates) within the operating conditions studied. The characteristics of drag reduction for different M. are similar, 
which shows it is not directly linked to L p , i.e., longer L p does not necessarily produce larger drag reduction. Figure 8 
depicts the unsteady nature of these LPM jets. The overall impact of the counterflowing jets on the flowfield near the nose 
of the body is itself shown in Fig. 9, through a comparison of the computed density gradient contours for the unmodified 
cone-cylinder and that of a truncated cone-cylinder which includes a counterflowing jet (M. = 4.02 ;P t j = 35.0 atm). 
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Figure 8. Mach number contours at different time instances for the cone-cylinder (Af. = 4.02). The x- andy- 
coordinates are given in inches. 
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Density Gradient Magnitude: 0 0.4 0.8 1.2 1.6 2 



Figure 9. Comparison of the computed density gradient contours for the unmodified cone-cylinder (top) and the 
truncated cone-cylinder, which includes a counterflowing jet (bottom). 


b. Effect of Jet/Nozzle Exit Diameter 

To help understand the effect of the jet/nozzle exit diameter (J.) on the jet penetration length, the nozzle with M. = 
2.94 was examined for three different nozzle exit diameters (d. = 1.27, 2.54, and 5.08 mm, or 0.05, 0.1 and 0.2 in., 
respectively) with the jet-to-base diameter ratio remaining constant (d. Id m = 1/4). The nozzle contour (bell-shape) was 
scaled to achieve different exit diameters. Given that the same nozzle shape was used in this study, the range of jet total 
pressures studied was roughly the same (3—15 atm) and the scaling accounted for the different mass flow rates. As shown 
in Figs. 10-12, L p increases with d. for a given d. Id m and M, but their ratio ( E p ld. ) is less sensitive to the change in d.. Also, 
for a given d., L p increases with NPR and jet flow rate to an optimal value and then decreases because the jet reverts back to 
SPM. The optimal NPR and jet flow rate condition for obtaining maximum L p varies slightly for the different nozzle exit 
diameters. Similar to the jet Mach number effect, the correlation between L Id. and NPR can be easily observed, compared 
to looking at the variation in L p with NPR alone. In addition, drag reduction increases with NPRs (and jet flow rates) within 
the operating conditions studied. The characteristics of drag changes for different d. are similar and show that drag reduction 
is not directly linked to L p . Though Fig. 11 reveals that a larger jet diameter yields higher drag reduction for the same jet 
total pressure (i.e., NPR), the main driver is the jet flow rate, as can be seen from Fig. 12. This effect occurs because a larger 
jet diameter corresponds to a higher jet flow rate for a given jet total pressure. Furthermore, longer L p in general produces a 
larger drag oscillation, which may be disadvantageous from a vehicle controllability perspective. 
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Figure 10. Effect of jet diameter on shock standoff distance with various jet total pressures for the cone-cylinder 
geometry (d. Id m = 1/4 and M. = 2.94). 



Figure 11. Effect of jet diameter on shock standoff distance and drag reduction with various NPR for the cone- 
cylinder geometry (d./d m = 1/4 and M. = 2.94). 




Figure 12. Effect of jet diameter on shock standoff distance and drag reduction with various jet flow rates for the 
cone-cylinder geometry (d./d m = 1/4 and M. = 2.94). 

c. Effect of Jet-to-Base Diameter Ratio 

Based on our previous study (Refs. 14 and 24), the jet-to-base diameter ratio ( d.ld m ) is the most significant factor 
affecting the formation of the LPM jet. The parameter d.ld m determines the height of the forward facing step above nozzle 
exit, which strongly influences the strength and pattern of flow recirculation in front of the truncated cone and the stability 


9 


of jet shear layer, and thus, the LPM jet formation. Three jet-to-base diameter ratios ( d . Id = 1/4, 1/8, and 1/16) were 
assessed with the same nozzle jet diameter d. = 0.635 mm (0.025 in) and jet exit Mach number (Af. = 2.94). The nozzle total 
pressure ranges from 5 to 15 atm, resulting in the same range of jet mass flow rate for different d.ld m due to the same d. and 
M.. The results indicate that, as the jet-to-base diameter ratio decreases, both the jet penetration length (L p ) and its ratio to 
the jet diameter (L p Id.) increase, as seen in Figs. 13-15. Furthermore, the effect of d.ld m on L p is larger than the effect of 
NPR (and jet flow rate), which confirms our observation from our earlier studies on the importance of this parameter. 24 
However, opposite to expectation, the drag increases as L p increases. The drag increase is believed to be caused by the 
additional friction loss induced by a larger base area associated with smaller d. Id m . This was confirmed by comparing the 
drag-reduction relative to the baseline and the corresponding truncated cone-cylinder geometries, shown in Fig. 16. When 
compared to the drag of the truncated cone without a jet, there is definite drag reduction in the presence of jet, implying that 
the jet does offset the additional drag induced by the front face of the truncated cone. The drag does reduce with d.ld m , when 
the drag reduction is computed in reference to the drag of the truncated cone without a jet. Hence, the overall drag increase 
or reduction depends on the magnitude of the pressure drag reduction due to the increase of L p ? and the magnitude of 
additional friction drag from the truncated face of the cone-cylinder. In this case, the addition of friction drag overwhelmed 
the pressure drag reduction, because the pressure drag is already very small for the given cone-cylinder geometry. Hence, 
for drag reduction, smaller d m (or larger d./d m ) should be used. The Mach contours in Fig. 15 once again demonstrate that 
the resulting counterflowing jets cannot regain the length lost in the truncation of the cone-cylinder, even with the longest 




Figure 13. Effect of jet-to-base diameter ratio on jet penetration length with various jet total pressures, for the cone- 
cylinder geometry (Af . = 2.94, and d. = 0.635 mm or 0.025 in). 




Figure 14. Effect of jet-to-base diameter ratio on jet penetration length and drag reduction with various NPR, for the 
cone-cylinder geometry (Af . = 2.94, and d. = 0.635 mm or 0.025 in). 
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Figure 15. Instantaneous Mach contours to show the effect of jet-to-base diameter ratio on jet development for the 
cone-cylinder geometry (Af. = 2.94, d. = 0.635 mm or 0.025 in). The x- and y -coordinates are given in inches. 




Figure 16. Effect of jet-to-base diameter ratio on drag reduction with respect to the baseline and the truncated cone- 
cylinders with various jet flow rates (Af. = 2.94, and d. = 0.635 mm or 0.025 in). 

To confirm our aforementioned observation that a smaller d m (or larger d.ld m ) should be used for drag reduction, two 
additional sets of cases were examined. The first set had a jet diameter of 1.27 mm (or 0.05 in.), and a base diameter of the 
same size, which provided the largest d.!d m (= 1). In the second set, the jet diameter was 0.635 mm (0.025 in.) and the base 
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diameter was 2.54 mm (0.1 in.), which gave d./d m = 1/4. The nozzle with an exit Mach number of 4.02 was used, and the 
jet total pressure ranged from 17 to 59 atm. As shown in Figs. 17-18, larger d.ld m produces shorter L , indicative of the 
stable SPM jet configuration, and a larger drag reduction. For smaller d. /d m , L p increases with P (and jet flow rate) to an 
optimal value and then decreases, which is a characteristic of LPM transition to SPM. Figure 19 reveals that at P { . = 35 atm, 
the counterflowing jet of d. .Id m = 1/4 is in LPM, while the jet structure of d. ,ld m = 1 exhibits the SPM characteristics. For the 
largest d.ld , drag reduction increases almost linearly with NPR and jet flow rate within the operating conditions studied. 

Thus, in the case of the cone-cylinder geometry, achieving LPM or longest L p is not necessary for optimal drag 
reduction. Furthermore, for the case of largest dj/d m , the penetration length of the counterflowing jet is still not long 
enough to compensate for the length lost in truncation. However, its shock standoff location is the one closest to the tip of 
the baseline cone-cylinder, as shown in Fig. 20, as the least amount of truncation in length was required to obtain this 
configuration. 




Figure 17. Effect of jet-to-base diameter ratio on jet penetration length with various jet total pressures, for the cone- 
cylinder geometry (Af . = 4.02). 



rates, for the cone-cylinder geometry (Af . = 4.02). 
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Figure 19. Instantaneous Mach contours showing the effect of jet-to-base diameter ratio on the jet development for 
the cone-cylinder (Af . = 4.02, P = 35 atm). 
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Figure 20. Comparison of Mach contours for the cone-cylinder for two different jet total pressures (Af . = 4.02, d.!d m = 
1). The x- and y -coordinates are given in inches. The nose tip of the unmodified cone-cylinder geometry was at x = 0. 


B. Quartic Geometry 

Given the sharp curvature near the nose of the quartic geometry, the largest nozzle that can be accommodated, 
without significantly altering the baseline geometry, corresponds to a nozzle with exit diameter of 0.7166 mm, or 0.0282 in. 
(d. Id m = 1/8). Computational results from the previous study 24 revealed that, establishing LPM jets with the d. Id m = 1/8 
nozzle was difficult, while there were no such difficulties with the d.ld = 1/16 nozzle. Based on these findings, the d.ld - 
1/8 nozzles were replaced with the d.ld m = 1/12 nozzles in this extended parametric study for the quartic geometry. Two 
nozzles with exit Mach numbers M. = 2.94 and 4.02, along with various jet total pressures were investigated. For M. = 4.02, 
the jet total pressure ranges from 17 to 47 atm, while the total pressure range is between 3 and 15 atm for the M = 2.94 case. 
All of the counterflowing jet configurations studied is summarized in Table 2. For the baseline quartic geometry (without a 
jet), the total drag, needed to estimate the percentage of drag change in the presence of counterflowing jet, was computed to 
be 0.48 N. 
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Table 2. Summary of counterflowing jet configurations studied for the quartic geometry 


d. 

j 

d.ld 

J m 

M 

j 

NPR 

P t,\ (atm) 

0.358 mm (0.0141 in) 

1/16 

2.94 

20.5 - 102.6 

3- 15 

0.358 mm (0.0141 in) 

1/16 

4.02 

116.3-321.5 

17-47 

0.478 mm (0.0188 in) 

1/12 

2.94 

20.5 - 102.6 

3- 15 

0.478 mm (0.0188 in) 

1/12 

4.02 

116.3-321.5 

17-47 
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Figure 21. Effect of jet-to-base diameter ratio on jet penetration length and drag reduction with various jet-to- 
freestream pressure ratios and jet flow rates for quartic geometry (Af. = 4.02, and d m = 5.733 mm or 0.2257 in). 


a. Effect of Jet-to-Base Diameter Ratio 

To assess the jet-to-base diameter ratio effect, we begin by comparing the results of d./d m =1/12 and 1/16 for M = 
4.02. As indicated by Fig. 21, a larger jet diameter (i.e., larger d. Id m ) leads to longer jet penetration length and less drag. 
However, the ratio of jet penetration length to jet diameter (L p Id/) for different d. is comparable. Similar to the cone- 
cylinder geometry, the use of the counterflowing jet on the quartic geometry exhibits the characteristics that L p increases 
over a range of NPR (and jet flow rate) to an optimal value and then decreases when the jets switch from LPM to SPM. The 
results of M. = 2.94 (Fig. 22) were similar in trend to those of M. = 4.02. The percentage of drag reduction for the quartic 
geometry was observed to be generally larger than that for the cone -cylinder geometry. This is because the implementation 
of a jet nozzle in the quartic geometry does not require any modification, thus, does not introduce any additional drag — 
resulting from increased normal surface area due to truncation. 
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Figure 22. Effect of jet-to-base diameter ratio on jet penetration length and drag reduction with various jet-to- 
freestream pressure ratios and jet flow rates for quartic geometry (Af. = 2.94, and d m = 5.733 mm or 0.2257 in). 


b. Effect of Jet Mach Number 

Two jet Mach numbers (M. = 4.02 and 2.94) were compared for d.ld m = 1/12 and 1/16, respectively, the results of 
which are summarized in Figs. 23-26. The results indicate that nozzles with higher jet exit Mach numbers can help establish 
jets with longer L p \ however, there was not a significant difference in the percentage of drag reduction. A higher jet Mach 
number is accompanied by a larger momentum loss in the nozzle. Hence, the reduction of drag due to longer shock stand- 
off distance is likely canceled out by the higher loss in the nozzle. In addition, the jet conditions (both total pressure and 
flow rate) to achieve the optimal L p vary for different jet Mach numbers, as shown in Figs. 23-26. 
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Figure 23. Effect of jet exit Mach number on jet penetration length and drag reduction with various NPR for quartic 
geometry ( d.!d m = 1/12, and d m = 5.733 mm or 0.2257 in). 



Figure 24. Effect of jet exit Mach number on jet penetration length and drag reduction with various jet mass flow 
rates for quartic geometry (d. Id m = 1/12, and d m = 5.733 mm or 0.2257 in). 



Figure 25. Effect of jet exit Mach number on jet penetration length and drag reduction with various jet-to- 
freestream pressure ratios and jet mass flow rates for quartic cone ( d.ld m = 1/16, and d m = 5.733 mm or 0.2257 in). 
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Figure 26. Effect of jet exit Mach number on jet penetration length and drag reduction with various jet-to- 
freestream pressure ratios and jet mass flow rates for quartic cone ( d.ld m = 1/16, and d m = 5.733 mm or 0.2257 in). 


V. Impact on Farfield Pressure Disturbance Signature 

Given that one of the main goals of this numerical study is to evaluate the effect of the counterflowing jets on the sonic- 
boom signature, we begin by examining the accuracy of the employed CFD solver in accurately predicting the farfield 
pressure signature of the baseline geometries without the counterflowing jet. This was done by comparing data from 
computational studies against those from experiments. 25, 26 Previous computational studies, 30 for predicting the sonic-boom 
signature of these two simple axisymmetric geometries have all been inviscid flow calculations, and have emphasized the 
need for specialized mesh generation or mesh-adaptation techniques to accurately capture the propagation of the pressure 
signature away from the object studied. However, capturing the farfield pressure signature in the presence of the 
counterflowing jets comes with severe constraints on the mesh requirements, for two reasons: (i) viscous effects are very 
important for capturing the development of the counterflowing jet and its associated shear layer, requiring a fine mesh near 
the wall and in the jet-development region; and (ii) the dimensions of the jet nozzle is very small compared to the 
computational domain size dictated by the distances at which experimental data were collected, resulting in a large mesh 
element count. Additionally, the employed flow solver does not have built-in solution-based mesh adaptation capability, 
which would have made the task of mesh generation easier. As a consequence, to identify the proper mesh required for our 
flow solver, Euler flow computations were performed for the baseline geometry of the cone-cylinder and the quartic 
geometry. Once, satisfactory results that compared well with available experimental data were obtained, a new mesh 
combining the farfield mesh requirement with the near-field mesh requirement, provided in Section III, was constructed. 

A. Baseline Results and Comparison Against Experiments 

The mesh for the baseline cases were generated by manually packing grid points around the Mach lines emanating from 
the nose and the expansion corner of the body, and iterating it a few times based on the solution obtained on each mesh. 
Since, the eventual goal was to study the influence of the jet on the farfield signature, the use of stretched, high-aspect ratio 
triangular cells along the shock lines was avoided, as that introduces excessive numerical dissipation in computations by the 
CESE based solver. For the cone-cylinder, the flow conditions (M x = 1.68, F too = 0.6215 atm, T too = 327.57 K) from Ref. 25 
were used, while the freestream conditions for the quartic geometry (M w = 1.41, P too = 0.6804 atm, T {ao = 311.0 K) match 
those of the test conditions in Ref. 26. 

The farfield pressure signature for the cone-cylinder geometry was calculated and compared at h/L =10 with different 
levels of grid refinement. The final mesh that provided good comparison against experimental data contained approximately 
4 million cells. The normalized pressure perturbation contour obtained from the use of this final mesh is shown in Fig. 
27(a). The comparison of the normalized pressure perturbation magnitude at hIL =10, obtained with the final mesh, against 
experimental data 25 is shown in Fig. 27(b). As observed from Fig. 27(b), the results from the final mesh are able to closely 
match the experimental data in the leading portion of the shock. In the expansion region, the computational results predicts a 
higher pressure (-13%) at the location of the lowest pressure, pointing to a need for further improvement of the mesh 
resolution in that region. 
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Figure 27. Baseline pressure disturbance (normalized) signature for cone-cylinder geometry. 

For the quartic geometry, the farfield pressure distribution was compared against that of the experiment at h/L - 5. 
The meshes were generated and refined in a manner similar to how it was done for the cone -cylinder geometry. The final 
mesh count was approximately 2.5 million cells. The normalized pressure perturbation contour and the comparison of its 
magnitude against that of the experiment at h/L = 5 are shown in Figs. 28(a) and (b), respectively. The source of oscillation 
of the wind tunnel measurements in the expansion region, -0.75 < x/L < 0.0 is unknown. The CFD results closely match the 
experimental data, in its prediction of the bow shock location and the angle of the expansion from the maximum 
overpressure. However, the CFD result once again predicts a higher pressure (-25%) at the location of the lowest pressure 
in the expansion region, similar to the outcome of the cone-cylinder geometry. 
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Figure 28. Baseline pressure disturbance (normalized) signature of quartic geometry. 

B. Impact of Counterflowing Jet on Farfield Pressure Signature 

a. Cone-cylinder Geometry 

To assess the impact of the counterflowing jet on the farfield pressure signature, the baseline pressure signature of 
cone-cylinder geometry without any geometry truncation or jet was computed once again, but under the freestream 
conditions of Mach 1.6 flow at an altitude of 13,700 m, because of its interest to the NASA’s FAP High Speed project. As a 
result, the mesh that was previously generated needed slight modifications, to accommodate the difference in the flow Mach 
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angle. A new mesh that included the appropriate mesh resolution in the jet development region and near -field region, along 
with adequate resolution in the farfield region (based on the baseline study) was constructed to study the impact of the jet on 
the farfield pressure signature. However, if the mesh distribution and computational domain were to be the same as those for 
the baseline study ( h/L =10), the resulting mesh would be so large that its use for a numerical parametric study would 
require computing resources and time exceeding the period of performance of this project. As a result, a reduced 
computational domain that allowed for comparison of pressure signal only up to h/L = 3 was used. 

The data for the baseline geometry and that with the counterflowing jet (M. = 2.94, P =9 atm., d. = 0.635 mm, 
and d. Id m = 1/4) was computed. The normalized pressure disturbance is shown at three different distances from the 
centerline (h/L = 1, 2, 3) in Fig. 29. As can be seen from Fig. 29, the jet does introduce additional jumps in the leading 
portion of the pressure rise, which can result in a favorable ground boom-signature with finite rise-time. However, it does 
not reduce the magnitude of the maximum pressure jump in the leading portion of the signature, in contrast to what was 
observed in the case of a blunt-body under higher-speed freestream flow conditions. 24 In addition, there is a significant 
impact on the expansion region of the signature, which could be a result of the inherent unsteadiness in the flow 
downstream of the body caused by the presence of the jet. To resolve the unsteadiness in the flow expansion region of the 
flow more accurately and understand the impact of counterflowing LPM jets on the pressure signature may warrant 
unsteady mesh-adaptation capabilities in the employed solver. The reason for the ineffectiveness of the counterflowing jet in 
reducing the farfield pressure signature for the cone-cylinder could be explained as follows. In the case of the baseline cone- 
cylinder, the leading shock is a weak oblique shock, which propagates at the Mach angle (determined by the freestream 
Mach number) well into the farfield. With the introduction of a counterflowing jet, the original cone geometry gets 
truncated, resulting in a leading bow-shock, which can be substantially weakened by the counterflowing jet. However, in the 
farfield region, the disturbance from the weakened leading shock, merges with a detached shock arising slightly 
downstream of the nose — a result of interaction between the jet recirculation flow and the expansion wave — thereby gaining 
in strength and propagating once again along the Mach angle. Moreover, with the length of jet penetration never being able 
to compensate for the length lost from truncation (i.e., the shock stand-off distance is smaller than the length of the cone 
truncated), delaying of the merging of those two shock systems becomes difficult. Overall, these results from the cone- 
cylinder are in line with those from Fomin et al., 23 who had used a similar type of geometry and concluded that 
counterflowing jets did not yield significant boom signature reduction for the cone-cylinder type geometry. 


Mach 2.94 Jet Jet_pt9_y25 (h/L=3) 

Jet total pressure = 9 atm Jet_pt9_yl8 (h/L=2) 



Figure 29. Comparison of pressure disturbance signatures at h/L = 1,2 and 3 with jet for the truncated cone-cylinder 
and and without the jet for the unmodified cone-cylinder geometry. 
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a. Quartic Geometry 

For the quartic geometry, the counterflowing jet condition (Afj = 4.02, P tj = 41 atm, dj = 0.358 mm, and dj ld m = 1/16) 
that produced the longest jet penetration length was used to assess the impact on farfield pressure signature. As seen from 
Fig. 30, the impact of the jet on the farfield pressure signature is shown at three different distances away from the centerline 
of the body ( h/L = 1, 2 and 3). Here, the presence of the jet does modify the sharp jump in pressure rise in the leading 
portion of the pressure signature, through the introduction of a more gradual increase in pressure jump. Such a modification 
will result in a favorable ground boom- signature profile with finite rise-time. As expected, in the near-field region {h/L = 1), 
the counterflowing jet has considerable impact on the pressure signature (reduction in pressure spike magnitude). However, 
as we move further away from the body, the degree of impact weakens; but the signature continues to display a favorable 
profile. However, whether or not the reduction in the magnitudes of pressure spike seen at /z/L= 2 and 3 continues further 
away from the body, needs further investigation. The investigation could be carried out by considering a larger 
computational domain, with additional solution-based mesh refinement or through the use of pressure signature propagation 
codes. 

Based on the results shown in Figs. 29 and 30, the counterflowing jet does not appear to be effective in reducing the 
farfield pressure signature when operating for low supersonic slender body configurations. This is in contrast to its behavior 
under moderate-to-high supersonic freestream conditions in the case of blunt -body geometries, where a strong bow shock is 
present. In the case of low supersonic slender body configurations, the initial shock system is already a weak one and hence, 
its alteration by the counterflowing jet is not sufficient to attenuate the pressure disturbance in the farfield. 
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Figure 30. Comparison of pressure disturbance signatures at h/L = 1,2 and 3 with and without jet for the quartic 
geometry. 


VI. Summary 

An in-depth parametric study, involving numerical computations, was carried out for investigating the effect of 
various counterflowing jet configurations on two supersonic slender-body models (cone-cylinder and quartic body of 
revolution). The primary objectives of the study were to obtain a better understanding of the relationship between the 
counterflowing jet’s penetration length and reduction of drag and boom signature for these two supersonic slender-body 
configurations. Different jet flow rates, Mach numbers, nozzle jet exit diameters and jet-to-base diameter ratios were 
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examined as part of this investigation, resulting in more than 100 cases that were computed. The results obtained from the 
parametric study, concerning conditions favorable to obtaining LPM jets, prolonging jet penetration lengths and their impact 
on drag reduction, can be summarized as: (1) the jet penetration length increases with NPR before it reaches an optimal 
value, beyond which the jet switches back from LPM to SPM; (2) though additional pressure and momentum drag is 
introduced by the implementation of nozzle jet, its magnitude is surpassed by the reduction of drag due to weakening of the 
shock by the impinging counterflowing jet for all the cases with the quartic geometry and most of the cases with the cone- 
cylinder; (3) increase of dj leads to longer L p , smaller drag force and larger percentage of drag reduction; (4) increase of 
NPR and jet flow rate leads to decrease of drag and increase of drag reduction percentage; (5) for the cone-cylinder 
configuration, achieving LPM or longest L p does not necessarily produce optimal drag reduction; (6) to achieve the longest 
L p , one should employ the smallest djld m , highest Mj, and the optimal NPR to yield LPM; (7) to obtain larger drag reduction, 
one should go towards a larger dj ld m and higher NPR with proper Mj\ and (8) longer L p implies a larger oscillation in drag, 
making it disadvantageous from a vehicle stability perspective. For drag reduction, the use of a counterflowing jet on the 
quartic geometry is more effective than on a cone-cylinder. For the cone-cylinder, not only is the counterflowing jet unable 
to compensate for the length of the cone lost from its truncation to facilitate the jet nozzle implementation, but the additional 
drag induced by the surface normal to the freestream present from the cone truncation, also decreases (for some cases it may 
even reverse) the effect of drag reduction benefit from the counterflowing jet. Whereas, for the quartic geometry, the 
counterflowing jet provides consistent drag reduction, for the various jet conditions and nozzle configurations studied. This 
is mainly because the implementation of a jet nozzle in the quartic geometry does not increase the surface area exposed to 
the freestream and thus, no additional drag is induced. 

Based on the analysis performed to study the impact of the counterflowing jet for the farfield pressure signature of 
the cone-cylinder and quartic geometry, in the case of a cone-cylinder geometry, even the longest jet penetration length did 
not result in a subsequent drop in the pressure disturbance levels continuously to the edge of the computational domain. For 
the quartic geometry, the pressure signature does drop up to a distance of three times the body length away from the 
centerline of the body. However, the significance of the reduction in pressure disturbance signature tends to decrease further 
away from the body, therefore, more investigations are needed to quantify the pressure spike magnitude in the farfield. For 
both geometries, the counterflowing jet does introduce a gradual jump in pressure in the leading portion of the pressure 
signature, which can lead to a ground signature with finite rise-time for the leading portion of the “N-wave” signal — 
beneficial to sonic-boom mitigation. One of the primary reasons for the counterflowing jet not being able to substantially 
reduce the pressure disturbance signature could be due to the low supersonic freestream conditions and the slender 
configurations studied. These conditions result in only a weak shock system that is not significantly impacted by the 
penetrating counterflowing jets, as opposed to what was seen in the case of blunt -body geometries operating under higher- 
speed freestream conditions. The result of farfield pressure signature comparisons also indicates the need for further 
improvement in computations. Currently, some activities are being performed in improving the accuracy of numerical 
computations with the CESE solver, such as the development of a high-order scheme and the use of solution-based mesh 
adaptation to achieve desired accuracy without using extremely fine mesh for the entire computational domain, which can 
lead to an undesirable computational cost. Hence, the effects of the counterflowing jet on the boom signature need further 
investigation in the future with the use of a higher-order CESE method and/or solution-based mesh adaptation. Based on the 
preliminary analysis of the impact of the counterflowing jet on the farfield pressure signature and more detailed analysis on 
the drag-reduction potential, this concept could still have aerodynamic benefits for quartic geometry type (blunt nose) 
objects in supersonic flight conditions, especially during the vehicle’s acceleration/deceleration phases. 
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